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ABSTRACT 

The opacity due to photodissociation of ^^MgH is investigated in the atmo- 
spheres of cool stars. The lowest two electronic transitions A^Il *~ X ^S"*" and 
B' ^S"*" ^ X are considered where the cross sections for the latter were 
published previously (Week, Stancil, & Kirby 2002a) while the former are pre- 
sented in this work. Model atmospheres calculated with the PHOENIX code are 
used to investigate the effect of the photodissociation opacity on spectra of cool 
stars. The A ^Tl ^ X ^S"*" photodissociation cross sections are obtained us- 
ing a combination of ah initio and experimentally derived potential curves and 
dipole transition moments. Partial cross sections have been evaluated over the 
accessible wavelength range A A 1770 — 4560 A for all rotational transitions from 
the vibrational levels v" = — 11. Assuming a Boltzmann distribution of the 
rovibrational levels of the X state, LTE photodissociation cross sections are 
presented for temperatures between 1000 and 5000 K. Shape resonances, aris- 
ing from rotational predissociation of quasi-bound levels of the A state near 
threshold, characterize the LTE photodissociation cross sections. A sum rule is 
proposed as a check on the accuracy of the photodissociation calculations. 

Subject headings: molecular processes — molecular data — stars: atmospheres 
— stars: late-type 
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1. Introduction 

The lack of accurate and complete molecular line and continuum opacity data has been 
a serious limitation to developing atmospheric models of cool stars (M and later), solar sys- 
tem planets, and Extrasolar Giant Planets (EGPs). Sophisticated modeling programs, such 
as PHOENIX (Hauschildt & Baron 1999), require high quality opacity data in order to pro- 
duce synthetic spectra and predict physical parameters (e.g., surface chemical composition, 
effective temperature, etc). Typically, atmosphere models include molecular bands with 
hundreds of millions of spectral lines, mostly derived from molecular band or Hamiltonian 
models. Moreover, atmospheric models do not consider the effect of molecular photodis- 
sociation processes, which may play a role in the opacity at visible and UV wavelengths. 
An exception is the work of Kurucz, van Dishoeck, & Tarafdar (1987) who included LTE 
photodissociation calculations of CH and OH in solar and cool stellar atmosphere models. 
Further, it has been suggested by Short & Lester (1994, 1996) that a missing opacity source 
between 2500 and 4000 A might be due to MgH, based on the ground rovibrational pho- 
todissociation cross section of Kirby et al. (1979) which peaks at 2920 A , but has a threshold 
at 3100 A. However, the LTE B' ^ X MgH calculations of Week et al. (2002a) 
find that the cross sections have significant amplitude to nearly 4500 A. 

In the present work, we have performed extensive calculations of rovibrationally-resolved 
photodissociation cross sections of ^^MgH through the A ^H X transition, using the 
most accurate available molecular data. Calculations were performed for the full range of 
313 rovibrational levels {v", J") in the ground electronic state. Assuming a Boltzmann dis- 
tribution of the rovibrational levels of the X state, LTE photodissociation cross sections 
are also presented for temperatures between 1000 and 5000 K. These photodissociation cross 
sections as well as those of Week et al. (2002a) for the B' ^ X transition have 
been included in cool stellar atmosphere calculations to test if MgH photodissociation is the 
missing opacity postulated by Short & Lester (1994, 1996). In this paper we investigate the 
importance of the bound-free opacity in both very cool dwarfs and the objects discussed in 
Short & Lester (1994, 1996) (K type giants and solar type dwarfs). 

Atomic units are used throughout to discuss the molecular calculations unless otherwise 
stated. 



2. Theory of Photodissociation 

Potential curves of the X and A ^H electronic states of MgH calculated by Saxon, 
Kirby & Liu (1978), as well as the connecting dipole transition moment, were used in the 
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present study. Shifts have been apphed to the potential curves, in order to correct for dis- 
crepancies between the theoretical dissociation energy and the experimental values given by 
Balfour & Cartwright (1976) and Balfour & Lindgren (1978). Details about this procedure, 
together with the fits used for short- and long-range internuclear separations, have been 
described in a previous publication (Week et al. 2002b). 

For absorption into the rovibrational continuum of electronic state / from the rovibra- 
tional state v" J" of electronic state i, the partial rotational photodissociation cross section 
can be written as (Week et al. 2002a) 

^ J"+i 

a^„j,. ^ 2.69 X 10-^' ^Ek'j>y>j» ^j,, ^ Sj'\Dk'j',v"j"\^ (1) 

j'=j"-i 

where AEk',j'^v"j" is the energy of the absorbed photon and Sj' are the Honl-London factors 
defined, for the yl X transition as 

r (J" - l)/2, J' = J" - 1 (P branch) 
Sj>{J") = I (2J" + l)/2, J' = J" (Q branch) (2) 
( {J" + 2) /2, J' = J" + 1 (R branch) 

according to the convention proposed by Whiting & NichoUs (1974). In equation (1), the 
term -D^Jj, ,y„j„ =< Xk'j'\D-^^\xv"j" > is the matrix element of the electric dipole transition 
moment responsible for absorption from the i into the / electronic state. The continuum 
wave functions Xk'j'{R) are normalized such that they behave asymptotically as 

Xfe'j'(i?)~sin(A:'i?-fJ' + 5jO. (3) 

The continuum and bound wave functions Xk'j' and Xv"j", respectively, were obtained 
by solution of the radial Schrodinger equation for the A and X states, respectively, using a 
standard numerical method (Cooley 1961). Wave functions were obtained on a grid with a 
stepsize of 1 x 10~^ a,o over the internuclear distance range 0.5 ^ R ^ 200 a,o- 

3. Photodissociation Calculations 

In Figures 1-3 are shown some representative rovibrationally-resolved cross sections for 
photodissociation through the A ^11 <— X transition^. 



^The complete list of data for A <— X photodissociation cross sections is available online at the UGA 
Molecular Opacity Project Database website http://www.physast.uga.edu/uganiop/ 
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Figure 1 presents our results for the partial cross sections a^/'j/', as a function of wave- 
length, for the transitions from the vibrational level v" = with the rotational quantum 
numbers J" = 0, 10, 20, 30 and 44 of the X ^E"*" electronic state. From the quantitative 
point of view, cross sections are seen to be small, if compared to the values obtained for pho- 
todissociation through the B' ^ X which are several orders of magnitude larger 
for the same wavelength range (Week et al. 2002a). This could be expected since the bound 
states of the A ^11 and X electronic states largely overlap, i.e., the sum of the Franck- 
Condon factors over all the upper vibrational levels v' is close to unity as originally pointed 
out by Kir by et al. (1979). Since the photodissociation cross sections are very small, a sum 
rule with rotational resolution has also been proposed as a check on the accuracy of our nu- 
merical calculations and is presented in the Appendix. On the other hand, some interesting 
qualitative features are observed, such as nearly periodic nodes and antinodes. This reflects 
the strongly oscillating behaviour of the continuum wave functions, since the A ^ X dipole 
transition moment decreases monotonically to zero as the internuclear separation increases. 
It is also worth noting the presence of shape resonances generated by predissociation of 
quasi-bound rotational levels near threshold. The magnitude of such resonances can reach 
several decades, as can be seen for J" — 44, the topmost rotational level for v" — 0. 

Partial cross sections for transitions from vibrational levels v" = 0, 4, 8 and 11 with the 
rotational quantum number J" = of the X electronic state are presented in Figure 2. 
As can be seen in the previous figure, the threshold hmit extends to larger wavelengths as v" 
and/or J" increase, corresponding to a decrease in the amount of energy necessary to reach 
the dissociation limit. In this connection, it is worth mentioning that the cross section value 
at threshold increases signiflcantly with the vibrational quantum number v". This may be 
explained by the fact that the broad maximum of the continuum wave function near the left 
turning point overlaps with the dipole transition moment over a wider range of internuclear 
distances. 

Figure 3 shows the partial cross sections for transitions from the topmost vibrational 
level v" = 11 with the rotational quantum numbers J" = 0, 1, 2 and 3 of the X 
electronic state. Few qualitative differences are observed between the maximum rotational 
number J" = 3 and J" = 0. In fact, in the case of such high-lying vibrational states, the 
spacing between all the rotational levels is small, thus resulting in few diferences between 
their respective wave functions which sample almost the same potential region. 

Assuming a Boltzmann population distribution of the rovibrational levels of the elec- 
tronic ground state X ^E"*", the LTE photodissociation cross section as a function of the 
wavelength is shown in Figure 4 for temperatures between T — 1000 and 5000 K. The 
LTE cross section shows strong shape resonances, arising from rotational predissociation 
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of quasi-bound levels near thresholds, as was observed previously for the partial cross sec- 
tions represented in Figure 1. The cross section exhibits two broad peaks centered around 
A = 4100 and 4400 A, respectively. Two broad peaks were also observed for the B' ^ X 
LTE cross sections, but centered near 2780 A and ~3800 A , with a mininum at 3200 
A (Week et al. 2002a). For the A ^ X transition, the relative minimum observed around 
A = 4200 A appears to be the result of a node common to most of the high-lying vibrational 
states, whose partial cross sections strongly contribute to the LTE photodissociation cross 
section. This feature is clearly illustrated in Figure 2 for the vibrational levels v" — 8 and 11 
on the wavelength range A A = 3800 — 4600 A, where the LTE cross section shows maxima. 
Let us note that the LTE cross sections increase with temperature, as expected, since more 
high-lying rovibrational states with lower thresholds are occupied at high temperature. 



3.1. Atmosphere Models 

The models used for this work were calculated as described in AUard et al. (2001). These 
models and their comparisons to earlier versions are the subject of a separate publication 
(AUard et al. 2001) and we thus do not repeat the detailed description of the models here. 
However, we will briefly summarize the major physical properties. The models arc based 
on the Ames H2O and TiO line lists by Partridge & Schwcnke (1997) and Schwenke (1998) 
and also include as a new addition the linelists for FeH by Phillips & Davis (1993) and for 
VO and CrH by R. Freedman (NASA- Ames, private communication). The models account 
for equilibrium formation of dust and condensates and include grain opacities for 40 species. 
In the following, the models will be referred to as "AMES-dusty" for models in which the 
dust particles stay in the layers in which they have formed and "AMES-cond" for models 
in which the dust particles have sunk below the atmosphere from the layers in which they 
originally formed. We stress that large uncertainties persist in the water opacities for parts 
of the temperature range of this work (Allard ct al. 2000). However, the MgH bound-free 
transitions are in the optical and only slightly affected by the quality of the water opacities. 

In addition to the opacities sources listed above and in Allard et al. (2001, and references 

therein) we added the new photodissociation opacity sources from this work and from Week 
et al. (2002a) to our opacity database. In order to assess the effects of the new MgH pho- 
todissociation data, we compare spectra calculated with and without these opacity sources. 
The original AMES grid was calculated for effective temperatures of M, L and T dwrafs. 
The hotter models in this work are based on the same physics as the AMES grid and merely 
differ in the effective temperature. The models used in the following discussion were all 
iterated to convergence for the parameters indicated. The high resolution spectra which 
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have the individual opacity sources selected are calculated on top of the models. The MgH 
photodissociation data turned out to be too weak to influence the temperature structure of 
the atmosphere. 

We calculated models with log(g')=5.0 and effective temperatures of 2000 K, 3000 K, 
and 4000 K as typical M and L dwarf parameters. To test if MgH is the missing opacity found 

by Short & Lester (1994, 1996), we calculated spectra with log(5')=4.5 and Tcg=5800 K (as 
a model to be close to a solar model) and with log(5') = 1.5 and Tcg=3900 (as a typical K 
giant). We found changes in all spectra of typically less then 2% when including the MgH 
photodissociation data. This is demonstrated in Figure 5 where the relative flux differences 
are shown. As can be seen, only the resonant dissociation transitions are strong enough to 
be visible. This is true for both the AMES-Dusty and AMES-Cond models although the 
absolute flux changes are smaller in the AMES-Dusty model since there is less flux in the 
optical due to strong dust opacity. The differences in the solar model are of the order of the 
numerical accuracy. We also tested a model with metallicity z = —1.0 (Teff=3000 K and 
log(5')=5.0) and found the effect due to the bound-free opacities to be practically unchanged 
compared to the models with solar abundances. 

In Figure 6 we show very high resolution synthetic spectra in the region of the strongest 
resonance with and without MgH bound-free data for the model with Tcfj=3000 K. The 
feature is due to B' ^ X photodissociation from the v" = 0, J" = 43 and v" = 1, J" = 40 
levels of the X state to quasi-bound levels of the B' state. However, the calculated wavelength 
of 4403 A is uncertain due to uncertainties in the long-range behavior of both the X and 
B' potential energy curves. 

4. Conclusion 

Photodissociation cross sections have been calculated for the A ^ X tran- 
sition of ^"^MgH for all rotational transitions from the vibrational levels v" = — 11 and 
over the accessible wavelength range A A 1770 — 4560 A. As predicted earlier by Kirby et 
al. (1979), using the Franck-Condon picture, photodissociation cross sections through the 
A are several orders of magnitude smaller than through the B' X pathway. 

However, interesting features such as a large number of shape resonances, arising from ro- 
tational predissociation near thresholds, appear in the LTE photodissociation cross sections 
calculated for temperatures between 1000 and 5000 K using a Boltzmann distribution of 
the rovibrational levels of the X state. Inclusion of the ^^MgH A ^ X and B' ^ X 
photodissociation cross sections in atmosphere models of cool stars results in only minor 
changes of the computed spectra. Though it does not explain the missing opacity found 
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by Short & Lester (1994, 1996), the MgH bound-free opacities will have to be considered 
for the analysis of high resolution spectroscopy at high S/N. Further, as MgH is a trace 
molecule, the opacity due to photodissociation may be more significant for molecules with 
larger abundances such as H2O. 
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National Science Foundation, and on the IBM SP of the NERSC with support from the 
DoE. We thank all these institutions for a generous allocation of computer time. This 
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Institute for Theoretical Atomic & Molecular Physics at Harvard University and Smithsonian 
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A. Appendix 

As a check on the accuracy of our numerical calculations, the sum rule for vibrational 
matrix elements proposed by Stephens & Dalgarno (1972) is extended to rotational resolu- 
tion. If {Xv'j',Xk'j'} is a complete set of discrete, bound rovibrational wave functions Xv'j' 
and energy normalized continuum wave functions Xk'J', and D an operator, then, from the 
equality case of the Schwarz inequality, we obtain 

< Xv"j"\D'\xv"j" >= ^J^PI'^P + T^iZ'i (Al) 
with, for transitions between bound states of the X and A ^H electronic states, 

VlTpf' = ^>S'j/| < Xv'j'\D\xv"j" > 1^ 

v',J' 

— < Xv'j"+i\D\Xv"j" > r 

v' 

+ -^j'l < Xv'j"\D\xv"j" > 1^ + S^,\ < Xv'j"-i\D\xv"j" > 1^} (A2) 



and, for transitions from bound states of the X electronic state to continuum states of 
the A electronic state, 



— / 'S'j'l < Xk> j'\D\xv" J" > \^dEhiji 

< Xk'j"+i\D\xv"j" > 1^ 
+ S^,\ < Xk'j"\D\xv"j" > \^ + Sj,\ < Xk'j"-i\D\xv"j" > \^}dEk', (A3) 



where we have introduced the Honl-London factors Sj,, Sj, and Sf, defined in equation (2) 
for the P-. Q- and i?-branches, respectively. 

For absorption from the rotational levels of the X ^E"*" {v" — 0) vibrational state, the 
sum rule, using the photodissociation cross section presented here and the bound-bound 
oscillator strengths of Week et al. (2002b), was satisfied to within an error of 0.0008% for 
J" = (D^/q = r'i'oV^^o™"'' = 8 X lO"*') up to a maximum error of 0.032% for the highest- 
lying rotational level, J" — 44 {V^^^ = 6 x 10~^). As the vibrational and rotational quantum 
numbers of the ground electronic state increase, the error increases, for example, with errors 
of 0.0078% for v" = 2 and J" = {v(^o = 8 x 10"^) and 7.08% for the topmost rotational 
level of v" = 2, J" = 39 {vl% = 7 x lO'^). 

As a remark, let us note that the present sum rule is aimed at checking the accuracy of 
the numerical calculations and in no way reflects the validity of the potential surfaces and 
the dipole transition moment. 
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Fig. 1. — Partial cross sections (Jy"j" for transitions from the ground vibrational level v" = 
with the rotational quantum numbers J" = 0, 10, 20, 30, and 44 of the X ^E"*" electronic 
state of ^"^MgH. On the left: linear scale; on the right: logarithmic scale. 
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Fig. 2. — Partial cross sections Gynj,, for tran- 
sitions from vibrational levels v" — 0, 4, 8, 
and 11 with the rotational quantum num- 
ber J" = of the X electronic state of 
24MgH. 
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Fig. 3. — Partial cross sections cr^" j" for tran- 
sitions from the vibrational level v" = 11 with 

the rotational quantum numbers J" = 0, 1, 2, 
and 3 of the X ^2+ electronic state of ^^MgH. 




Fig. 4. — The MgH LTE photodissociation cross section as a function of wavelength for 
temperatures between T — 1000 and 5000 K. On the left: linear scale; on the right: loga- 
rithmic scale. 
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Fig. 5. — Relative flux differences between model spectra after including MgH photodissoci- 
ation data. The parameters are indicated. All models are AMES-Cond, except the top left 
one. 
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Fig. 6. — Synthetic spectra at 0.001 A res- 
olution in the region of the bound-free res- 
onance. The models have Teff=3000 K and 
log(g)=5.0. The black line is without bound- 
free absorption, the grey line is with bound- 
free absorption. 



